How to make graphene superconducting 
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Graphene [l[ is the physical reahzation of many funda- 
mental concepts and phenomena in solid state-physics [2], 
but in the long list of graphene remarkable proper- 
ties a fundamental block is missing: supercon- 
ductivity. Making graphene superconducting is relevant 
as the easy manipulation of this material by nanolyto- 
graphic techniques paves the way to nanosquids, one- 
electron superconductor-quantum dot devices 0, S], su- 
perconducting transistors at the nano-scale[^ and cryo- 
genic solid-state coolers [lol|. 

Here we explore the doping of graphene by adatoms 
coverage. We show that the occurrence of superconduc- 
tivity depends on the adatom in analogy with graphite 
intercalated compounds (GICs). However, most surpris- 
ingly, and contrary to the GIG case[ll, 12|, Li covered 
graphene is superconducting at much higher temperature 
with respect to Ga covered graphene. 

As graphene itself is not superconducting, phonon- 
mediated superconductivity must be induced by an en- 
hancement of the electron-phonon coupling (A), 



A = 



N{0)D^ 



(1) 



In Eq. [T]iV(0) is the electronic density of states (DOS) at 
the Fermi level, D is the deformation potential, while M 
and ujph are effective atom mass and phonon frequency 
that in metallic alloys reflect the role of the different 
atomic species and phonon vibrations involved in su- 
perconductivity. In undoped graphene A is small and 
phonon-mediated superconductivity does not occur as 
the small number of carriers, intrinsic in a semimetal, 
leads to a vanishingly small N{0). In this respect the sit- 
uation is similar to the bulk graphite case, where, with- 
out intercalation of foreign atoms superconductivity is 
not stabilized. 

A first guess to induce superconductivity could then be 
to fill by rigid-band doping the carbon tt— states in order 
to have enough carriers. However, beside the fact that 
the TT— DOS grows very slowly with doping, this hurts 
against two major difficulties. First, even if the defor- 
mation potential related to coupling between tt— bands 
and in-plane phonon vibrations is very large and leads 
to Kohn anomalies these vibrations are highly en- 
ergetic {(jjph « 0.16 eV) and A is suppressed by the 
factor in the denominator. Second, symmetry forbids the 
coupling between tt— states and the softer out-of-planes 
vibrations. 

In order to promote coupling to carbon out-of-plane 
vibrations it is then necessary to promote new electronic 




FIG. 1. Crystal Structure of metal adatom covered 
graphene: yellow (blue) spheres represent carbon (metal) 
atoms. Adatoms sit on the hollow site of graphene layer. 



states at the Fermi level as it happens in GICs. Indeed 
in superconducting GIGs an intercalant bands (interlayer 
state) occurs at the Fermi level [14.] having multiple ben- 
eficial effects on A as (i) the number of carriers is en- 
hanced, (ii) coupling to G out-of-plane vibrations is pro- 
moted, (iii) coupling to intercalant vibrations occurs with 
a corresponding enhancement of the deformation poten- 
tial and a reduction of the effective Afcj^^ term in the 
denominator of Eq. [T] 

In GIGs not all kind of intercalant atoms leading to 
an interlayer state are equally effective in increasing Tc 
[lit ■ The larger Tc is indeed obtained when the distance 
(h) between the intercalant atom and the graphite plane 
is smaller. The reason is that the closer the intercalant 
electrons are to the planes, the larger are the carbon out- 
of-planes and the intercalant-modes deformation poten- 
tials D. This is demonstrated by the increase of Tc from 
BaGe (not superconducting) to SrCfi (T c = 1.65 K)[17| 
and finally to GaGg (Tc = 11.5 K)[ll|,[i3,[lll and by the 
increase of GaGg critical temperature under hydrostatic 
pressure [l3|. 

However a too small intercalant-graphite layer distance 
could also be detrimental for superconductivity. Indeed 
the quantum-confinement of the interlayer state in a too 
narrow region between the planes could result in a up- 
shift of the intercalant band well above the Fermi energy. 
In this case the ionization of the intercalant atom is com- 
plete and superconductivity is totally suppressed, as in 
bulk LiGeQ. 

Following these accepted guidelines for graphite inter- 
calated compounds, we explore the doping of graphene 
by foreign-atoms coverage (see Fig. [T]) and the possibility 
in in generating superconductivity in graphene by first- 
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TABLE I. Calculated physical parameters: Optimized struc- 
tural parameters (in A), electron-phonon coupling (A), loga- 
rithmic frequency average (cJjog in cm~^) and superconduct- 
ing critical temperature (T^ in Kelvin) for CaCe and LiCg 
systems. 
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h 


A 




Tc 


CaCe bulk 


4.30 


2.19 


0.68 


284.3 


11.5 


CaCs mono 


4.26 


2.24 0.40 


309.9 


1.4 


LiCe bulk 


4.29 


1.83 


0.33 


715.7 


0.9 


LiCe mono 


4.26 


1.83 


0.61 


277.8 


8.1 



principles density functional theory calculations lH l^l ■ 

The parallel with GICs is tempting, and in particu- 
lar, considering that Calcium intercalated graphite shows 
the highest superconducting critical temperature among 
GICs, the first example we consider is Calcium doped 
graphene and for comparison we simulated Ca interca- 
lated graphite. Calcium, as other alkaline metals, ad- 
sorbs in the hollow site of graphene (see Fig[Tl Tablc|l| 

The electronic band structure of both compounds is 
reported in Fig. [2] and shows that the deposition of Ca 
on-top of graphene still preserve the interlayer band at 
the Fermi level present in the bulk, fundamental for the 
appearance of superconducting phase in GICs. The re- 
moval of quantum confinement along the c-direction, due 
to the periodically stacked graphene layers, in passing 
from the bulk to the monolayer system, lowers the en- 
ergy of the IL state, that will result more occupied and 
carbon 7r-bands less doped. 

However, the sole presence of the IL state at the Fermi 
level being a necessary condition, is not sufficient to guar- 
antee large electron-phonon coupling, unless this state is 
not coupled with out-of-plane vibrations of graphene lat- 
tice. As the localization of the IL state as close as possible 
to the graphene layer 3 21, 22| enhances the coupling, 
in Figl3]we compare the planar-average of the IL charge 
density (calculated at the F-point of the Brillouin Zone) 
in the case of bulk and monolayer CaCg- 

In the monolayer, the IL charge density spills-out in 
the vacuum region while in the bulk case, it is much 
more confined in between the graphene and adatom lay- 
ers (roughly the region from z—Q A to z—2 A ). 

Based on this last consideration, we can infer that 
the superconducting critical temperature of CaCg mono- 
layer should be lower than its bulk counterpart. In 
order to verify this hypothesis, we calculated the vi- 
brational spectrum and the electron-phonon coupling 
A = 2/ duja'^F{ijj)/ijj of both compounds. The results 
are summarized in Table HI Fig|3]and FigjSl 

First of all we notice that the monolayer system is dy- 
namically stable, not showing the tendency to displacive 
instabilities, with a phonon dispersion characterized by 
three regions: the low energy region of adatom related 
modes (up to 300 cm~^) extending up to 400 cm~^ when 



> 1 

£? 

^-I 



-2 







\ Caq bulk 

















> 

>> 

a 
w 







\ CaC mono 

X \ — 





















> 1 

£? 
a 

W .1 



-2 

2 
1 







"V LiC, bulk 
















\ 




/ 


i 





> 



£? 
a 

W -1 
-2 







\ \\ LiC, mono 




\ 




/ 


1 



M 



K 



r 5 

DOS (states/eV) 



FIG. 2. Electronic band structure: From top to bottom: 
band structure of CaCe bulk compound (plotted along the 
high symmetry lines of the monolayer lattice), CaCs mono- 
layer and LiCe bulk and monolayer. Thickness of the bands 
is proportional to the intercalate/adatom s character. The 
Fermi energy is set at zero. 



mixed with carbon out-of-plane modes (C^), an interme- 
diate region of Cz modes (400-900 cm~^) and the high 
energy region characterized by C-C stretching modes. 

The major difference between bulk and monolayer 
cases is the softening of Ca vibrations. From the inspec- 
tion of the a'^F{uj) (which gives the contribution of each 
frequency to the total electron-phonon coupling), we note 
the low contribution of the Cz modes (around 500 cm~^ 
to the total electron-phonon coupling in the monolayer 
case with respect to the bulk CaCe, as expected from 
the above considerations. While the critical temperature 
of the bulk compound is 11.5 K (see Methods) in the 
monolayer case it lowers down to 1.4 K. 

This scenario suggests that, as a general rule, the re- 
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FIG. 3. Interlayer state wavefunction: Planar (in the xy 
direction ) average of |<^r,/Lp along the perpendicular (with 
respect to the graphene layer) direction (2). Vertical dashed 
lines represent z-position of Ca and Li adatoms. 



moval of quantum confinement should be detrimental for 
the electron-phonon coupling when compared with re- 
spect to same bulk (three dimensional periodic) system, 
due to the shift of the IL wavefunction away from the 
graphene layer. Every metal covered graphene should 
have (at least in the same stoichiometry) a reduced su- 
perconducting critical temperature with respect to the 
corresponding GIC. 

However, there is at least one case (to the best of 
our knowledge) among GICs that can be further ex- 
plored, namely stage- 1 Lithium intercalated compound, 
LiC6[23[, in which the IL state is completely empty (see 
Fig. [2]), as the strong confinement along the z direction 
(FiglS]) prevents its occupation and for this reason it is 
not superconducting. In this case, removal of confine- 
ment along c direction (Figl3]) should bring the IL at the 
Fermi level, as confirmed by the calculated band struc- 
ture (Figl2]). Inspection of Figl3] shows that the spatial 
extension of the IL for Li is the same as in Calcium, but 
being the IL strongly localized around the adatom and 
due its closer position to the graphene layer, we can ex- 
pect an enhancing of the total electron-phonon coupling. 

The phonon spectrum confirms, even in this case, the 
dynamical stability of the monolayer system and the com- 
parison with the bulk counterpart reveals that low en- 
ergy adatom modes and carbon vibrations in the direc- 
tion perpendicular to the plane are strongly softened in 
the monolayer. This suggests that the first two phonon 
branches, related to the in-plane displacements of Li 
atoms, and the undispersing Einstein mode displacing 
Li and C atoms out-of-phase along the z-direction at 
around 500 cm~^, should undergo enhanced electron- 
phonon coupling with respect to the bulk. On the con- 
trary, in the CaCg monolayer, adatom vibration and C 
vibrations along z are essentially at the same energy as 
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FIG. 4. Dynamical properties: Phonon frequencies disper- 
sion and phonon density of states (PDOS) of CaCe (upper 
panel) and LiCe (lower panel). Bulk (black) and monolayer 
(red) cases are considered resolving the adatom (shaded) con- 
tribution to the total (solid line) PDOS. 



in bulk CaCe. 

We found that bulk LiCe is a weak electron-phonon 
coupling A = 0.33 superconductor with an estimated su- 
perconducting critical temperature Tc = 0.9K^], theo- 
retically confirming the absence of superconductivity in 
LiCg. In the case of monolayer LiCg, the total electron- 
phonon coupling is A — 0.61, with a relevant enhance- 
ment of the superconducting critical temperature up to 
8.1 K. The main contributions to the electron-phonon 
coupling come from the low-energy Li modes and carbon 
vibrations along z, as expected from the above consider- 
ations, with a appreciable contribution from C-C stretch- 
ing modes (0.1). 

In summary, the ideal conditions for inducing super- 
conductivity in graphene are {i) bring the interlayer state 
at the Fermi energy (ii), localize it as close as possible to 
the graphene plane. This will ensure that the presence 
of the interlayer state effectively switch-on the dormant 
electron-phonon coupling of Cz modes that are inactive 
in the bulk. In addition, the in-plane displacement of 
adatoms will promote two additional scattering mecha- 
nisms, an intra-band contribution due to the interlayer 
state and an inter-band one due to the interlayer- tt scat- 
tering. Graphene can he made superconducting by the 
deposition of Li atoms on the top of it. 

As Li readily intercalates into graphite even at 100 
k[25} it is possible to incorporate Li atoms even below 
the graphene sheet. This possibility, based on our calcu- 
lations, is indeed favorable with respect to the developing 
of a superconducting phase. In fact, the double adsorp- 
tion should double the presence of the interlayer state at 
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FIG. 5. Eliashberg function: CaCs and LiCe (bulk and mono- 
layer) a'^F{uj). In the same graph is reported A(cj). 



Bulk LiCe has an aa stacking and was simulated with and 
electron-momentum mesh of 12x12x10 and 6x6x6 phonon- 
momentum grid for the calculation of phonon frequencies. A 
30x30x25 electron- momentum mesh was used for the electron- 
phonon coupling. 

The Eliashberg function a'^F{uj) is defined as: 



1 



N{0)Nk 



E 



ISkn, 



k:n,qm,i.' 

l5(£kn)'5{ek+qm)5('^ - l^q) (2) 

The total clcctron-phonon coupling X{u]) plotted in Fig[5]is defined 



A(a;) 



JO 



du)' 



(3) 



The total electron-phonon coupling is A(cj — ^ oo). The supercon- 
ducting critical temperature was estimated using the Allen-Dynes 
formula with /i* = 0.115, which fits the experimental critical tem- 
perature measured in CaCe GIC[ll||. 

Work supported by: CINECA-HPC ISCRA grant, by the EU 
DEISA-SUPERMAG project and by a HPC grant at CASPUR. 
Part of the calculation were performed at the IDRIS supercomput- 
ing center (project 91202). 



the Fermi level (one coming from each side of graphene) 
as we indeed verified. 

Li2C6 shows an increased coupling in all the frequency 
range due to the additional strongly coupled interlayer 
bands, with an electron-phonon coupling of 1.0 giving 
rise to critical temperature around 17-18 K. 

We conclude noting once more that bulk LiCe un- 
dergoes a remarkable metal-to-superconductor transition 
when exfoliated to one layer. The reason is that in a sin- 
gle layer LiCg the interlayer state cross the Fermi level, 
enhances the electron-phonon coupling and induces su- 
perconductivity. Our work demonstrates that supercon- 
ducting properties of adatoms on graphene are very dif- 
ferent from their bulk GIG counterparts. 

METHODS 

The results reported in the present paper were obtained from 
first-principles density functional theory in the local density 
approximation The QUANTUM-ESPRESSO O package was 
used with norm-conserving pseudopotentials and a plane-wave cut- 
off energy of 65 Ry. All the structures considered were relaxed to 
their minimum energy configuration following the internal forces 
on atoms and the stress tensor of the unit cell. 

The monolayer systems were simulated in the \/3 X \/3^?30° 
in-plane unit cell (with respect to standard graphene lattice, see 
FiglTJ with one adatom per unit cell. Phonons frequencies were 
calculated in the linear-response technique on a phonon wavevec- 
tor mesh of 12x 12 with a 14x 14 uniform electron momentum grid. 
The electron-phonon coupling parameter was calculated with elec- 
tron momentum fc-mesh up to 40x40 

CaCe bulk compound was simulated in the experimentally 
found structure with a/37 stacking [Till with a uniform electron- 
momentum fc-mesh integration of 8 X 8 X 8. The phonon frequen- 
cies were calculated on a 4 X 4 X 4 phonon-momentum mesh and 
the electron-phonon coupling integrated on a 20 X 20 X 20 electron- 
momentum mesh. 
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